Stress— softening and recovery of elastomers 

Aleksey D. Drozdov^ and Al Dorfmann^ 

^ Institute for Industrial Mathematics, 4 Hanachtom Street 

84311 Beersheha, Israel 
^ Institute of Structural Engineering, 82 Peter Jordan Street 
1190 Vienna, Austria 



Abstract 

A constitutive model is developed for the mechanical response of elastomers at fi- 
nite strains. A polymer is treated as a network of linear chains linked by permanent 
(chemical crosslinks) and temporary (entanglements and van der Waals forces) junc- 
tions. Temporary junctions are assumed to be in two states: loose (passive) when 
they impose only topological constrains on available configurations of chains, and 
tight (active) when their effect is tantamount to that for crosslinks. Stretching of 
a specimen implies that some loose junctions become active, which decreases the 
average length of a chain. 

A long chain is treated as an ensemble of inextensible strands connected in sequel. 
Two neighboring strands are bridged by a bond which may be in two conformations: 
flexed (trans) and extended (cis). A bond in the flexed conformation is modeled 
as a linear elastic solid, whereas the mechanical energy of a bond in the extended 
conformation (two rigid rods directed along a straight line) is disregarded. For a 
virgin specimen, all bonds are in the flexed conformation. Under loading some bonds 
are transformed from flexed to extended conformation. 

Stress-strain relations for a rubbery polymer and kinetic equations for the trans— 
cis transition are derived using the laws of thermodynamics. Governing equations 
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are determined by 5 adjustable parameters which are found by fitting experimental 
data in uniaxial tensile tests on natural rubber vulcanizates with various amounts of 
crosslinks. Fair agreement is demonstrated between results of numerical simulation 
and observations with the elongation ratio up to /c = 8. We analyze the effects of 
cyclic loading, thermal annealing and recovery by swelling on the material constants. 

Key— words: Elastomers, Stress-softening, Recovery, Rigid-rod chains. Temporary net- 
works 



1. Introduction 

This study is concerned with modeling stress-softening and subsequent recovery (by anneal- 
ing at an elevated temperature and by swelling) of elastomers. Stress-strain relations for 
unfilled and particle reinforced rubbery polymers have been the focus of attention in the 
past four decades. This may be explained by numerous applications of rubbery-like materi- 
als in industry (vehicle tires, seals, shock absorbers, flexible joints, etc.), as well as by some 
peculiarities in their mechanical behavior whose physical mechanisms are not quite clear. A 
renewal of the interest to the mechanical response of elastomers observed in the past decade, 
see Refs.llHll, is associated with the development of constitutive equations that (i) account 
for the molecular structure of polymers at the micro-level and (ii) describe some features of 
their response induced by time-dependent changes in the internal structure (viscoplasticity, 
damage, micro- fracture, etc.). 

Among these peculiarities, one of the most important for applications is the MuUins 
effect^. This phenomenon is evidenced in tensile0H^, compressive^ and shear0 tests as a 
noticeable difference between the stress-strain curve for a virgin specimen and that for the 
material reloaded after retraction. A conventional standpoint is that cyclic pre-stretching 
with 3 to 6 cycles and the maximum elongation ratio of about k = 2 makes the response 
of an elastomer repeatable§. However, the number of cycles, as well as their amplitude and 
frequency are chosen by the trial-and-error method and substantially depend on what the 
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repeatability of a stress-strain diagram means for the experimentalist. 

Although the physical mechanism of the Mullins effect in unfilled elastomers remains 
hitherto obscure, it is traditionally associated with material damage00'El, non-affinity of 
deformation for a polymeric network0'0, viscous drag excerted on the network chains by 
their environmentii and mechanically-induced crystallization of rubbery polymers^. These 
reasons may, however, be questioned, because observations reveal that with an increase in 
strains above the maximal level reached at pre-stretching, the stress-strain curve returns to 
that for a virgin sample. 

The present study aims to explain stress-softening of elastomers under cyclic loading 
by an increase in the average size of globules formed by polymeric chains (uncoiling of 
macromolecules) . For this purpose, we develop a phenomenological model for the stress- 
strain response of rubbery polymers and determine its parameters by fitting experimental 
data in uniaxial tensile tests for virgin, pre-stretched and recovered specimens. 

An elastomer is taken as a network of linear macromolecules bridged by permanent 
(chemical crosslinks) and temporary (physical crosslinks, entanglements and van der Waals 
forces) junctions. For the sake of simplicity, it is presumed that junctions move affinely with 
the bulk material. The novelty of our approach is that the number of junctions is not fixed, 
but assumed to alter under loading. To explain this phenomenon, we ascribe two possible 
states to an entanglement: loose (passive) and tight (active). An entanglement in the loose 
state does not create a junction, but imposes only topological constrains on the number of 
available configurations for a chain. On the contrary, an entanglement in the tight state is 
thought of as the same junction as that formed by a chemical crosslink. The transition of 
an entanglement from its passive state to the active state (tightening of entanglements) is 
supposed to be driven by mechanical factors. Introducing a simple kinetic equation for the 
rate of this transition, we study the response of a temporary network with a time-varying 
number of junctions. 

A chain that connects two neighboring junctions is modeled as a sequence of strands 
(statistically independent segments) bridged by bonds^. With reference to the Kratky- 
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Porod conceptEj, strands are taken as inextensible (rigid) rods, whereas the mechanical 
energy of a chain is determined as the sum of strain energies for bonds. It is postulated that 
a bond may be in two stable conformations: flexed (trans) and extended (cis). A bond in the 
flexed conformation is thought of as a linear elastic solid, whereas the mechanical energy of 
a bond in the extended conformation vanishes. Constitutive equations for an elastomer and 
the kinetic equations for the transition of bonds from their flexed to extended conformations 
are derived by using the laws of thermodynamics. The development of stress-strain relations 
is based on the hypothesis that the characteristic time for changes in the concentration of 
active entanglements substantially exceeds that for mechanical deformation. 

The paper is organized as follows. Deformation of a chain is discussed in Section 2. 
The mechanical energy of a temporary network is determined in Section 3. Stress-strain 
relations are developed in Section 4 using the laws of thermodynamics. Kinetic equations 
for transition of entanglements from their loose to tight state are introduced in Section 5. 
Uniaxial extension of a specimen is studied in Section 6. Results of numerical simulation 
are compared with experimental data for natural rubber vulcanizates in Section 7. Section 
8 deals with the effects of cyclic stretching and recovery on experimental constants. Some 
concluding remarks are formulated in Section 9. 

2. Deformation of a long chain 

An elastomer is modeled as an ensemble of linear macromolecules bridged by crosslinks, en- 
tanglements and van der Waals forces. Active links between macromolecules are thought of 
as junctions (permanent in the case of chemical crosslinks and temporary for entanglements 
and van der Waals forces). A sequence of mers (belonging to a polymeric molecule) between 
two subsequent junctions is associated with a chain. Chains are divided into strands (statis- 
tically independent segments) bridged by bonds. In accord with the Kratky-Porod model, 
a chain is treated as an aggregate consisting of + 1 identical inextensible strands linked 
in sequel. The average number of bonds in a chain change under stretching of a specimen 
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because of slippage of entanglements with respect to chains and transition of entanglements 
from their loose to the tight state. This implies that is taken as a function of time, 
= N{t). 

We adopt a zigzag mode@ and assume that a bond bridging two neighboring strands is 
characterized by one of the two stable conformations: flexed and extended. In the stress- 
free state the angle between two strands linked by a bond in the flexed conformation equals 
6 G (0, tt) (the same for all bonds) and that for a bond in the extended conformation equals 
TT (see a sketch depicted in Figure 1 of Ref.0). In thermal equilibrium before deformation 
all bonds are in the flexed conformation. Mechanical loading activates chains, which implies 
that conformations of some bonds alter. The numbers of bonds in flexed and extended 
conformations, iVf(t) and N^it), obey the balance law 

Ni{t) + N,{t) = N{t). 

Introducing the average concentration of bonds in the extended conformation, n = n(t), we 
obtain 

Ni{t) = N{t)[l - n{t)], iVe(t) = N{t)n{t). (1) 

In a deformed state the angle between strands linked by a bond in the flexed conforma- 
tion changes, whereas for a bond in the extended conformation (modeled as two rigid rods 
directed along a straight line) this angle remains unaltered. The strain, e, from the stress- 
free state of a chain to its deformed state equals the sum of strains for bonds in the flexed 
conformation, 

e = A^fCf. 

This equality together with Eq. (1) implies that 

Bonds in the flexed conformation are modeled as linear elastic solids with the mechanical 
energy 
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whereas the strain energy of bonds in the extended conformation vanishes. The mechanical 
energy of a chain, w, equals the sum of the mechanical energies for individual bonds. It 
follows from Eqs. (1) and (2) that 
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^^^^ = 2N^y 

To express the strain in a chain, e, in terms of the macro-strain tensor for a network, e, we 
consider a chain connecting neighboring junctions Ai{t) and A2{t). Let foi(t) and fo2(t) be 
radius vectors of these points in the stress-free state and fi{t) and f2{t) their radius vectors 
in the deformed state at time t > 0. The end-to-end length of the chain reads So{t) in 
the stress- free state and S{t) in the deformed state. Introducing the guiding vector in the 
stress-free state, I, (the unit vector directed along the end-to-end vector for the chain), we 
obtain 

ro2{t) - foi{t) = So{t)l. 

In the deformed state, the junctions occupy points with the radius vectors 

fl{t) = foi(t) + U{t, foi(t)), f2{t) = fo2(t) + U{t, fo2(t)), 

where u(t,f) is the displacement vector at point f for transition from the stress- free state 
to the deformed state at time t. The end-to-end vector for the chain in the deformed state 
is given by 

R{t) = f2{t) - n{t) = So{t)l + [u{t,foi{t) + 5o{t)l) - u{t,foi{t))]. 

Neglecting terms beyond the first order of smallness compared to 5o(^), we find that 

R{t) = 5o{t)l- [/ + VoM(t)], 

where Vq is the gradient operator in the stress-free state, / is the unit tensor, the dot stands 
for inner product and the argument foi is omitted. In terms of the radius vector in the 
deformed state, f, this equality reads 
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R{t) = So{t)l-Wof{t) = So{t)[Vof{t)r ■ I (4) 

where T stands for transpose. The end-to-end length of the chain in the deformed state, 
5{t), is given by 

5'' ^ R-R. 

This equahty together with Eq. (4) yields 

5\t;i)^5l{t)l-C{t)-l (5) 
where 

c{t) = Vof(i) ■ [v^mY (6) 

is the Cauchy deformation tensor for transition from the stress-free state of the network to 
its deformed state at time t. 

The extension ratio for the chain, X{tJ), is defined as the ratio of the current end-to- 
end length, S{t,l), to that of the chain in its activated stress-free state, S°{t,l). The latter 
state is defined as a state in which the chain is unloaded, but the numbers of bonds in 
various conformations coincide with their current values, iVf(t, J) and Ne{t, I). It differs from 
the equilibrium stress-free state of the chain, where the number of bonds in the extended 
conformation vanishes. The difference between the end-to-end lengths of a chain in the acti- 
vated state, 6°{t, I), and in the equilibrium state, So{t), determines the end-to-end elongation 
driven by transformation of bonds from their flexed conformation to the extended confor- 
mation. We assume that the transformation-induced end-to-end elongation is proportional 
to the number of bonds acquiring the extended conformation, 

S%tJ) = Som+Vn{t,l)], (7) 

where rj > characterizes an increment of the end-to-end length driven by an individual 
transition. The quantity rj reflects the average size of a globule created by a polymeric chain: 
when T) is small, the chain is rolled into a tight coil, whereas an increase in r] is tantamount 



to unfolding the coil. Because the parameter t] describes the current state of a chain, it 
should be thought of as a function of time (strain). To simplify calculations, we treat it as 
a constant, but suppose that rj can change from one stress-strain curve to another in cyclic 
tests. 

It follows from Eqs. (5) and (7) that the Hencky strain for a chain, e = In A, reads 
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/ • Cit) ■ I 



In 



l + rin{t,l) 



(8) 



5°(t,0 

Equations (3) and (8) determine the mechanical energy of a chain, w, in terms of the Cauchy 



deformation tensor for the network, C. 



3. Strain energy density of a network 

We adopt the conventional hypothesis that the excluded-volume effect and other multi-chain 
effects are screened for an individual chain by surrounding macromoleculesil. This implies 
that the energy of interaction between chains is neglected (under the incompressibility con- 
dition for the network) and the mechanical energy of the network equals the sum of the 
strain energies for individual chains. Assuming the distribution of chains to be isotropic, we 
find the concentration of chains in a network (per unit mass) with guiding vector /, 

X(t,l) = ^ sin MMif, 

4:71 

where and cp are Euler's angles which determine the position of the unit vector / with 
respect to some Cartesian coordinate frame and H(t) is the number of chains (per unit 
mass) at time t > 0. To determine the strain energy of a network, we multiply the number 
of chains, X{t,l), by their mechanical energy, w, and sum the results for various guiding 
vectors, I, 

^' 87f N{t) Jo ^ Jo l-n(t,t9,(^) ^' 
The quantities N(t) and S(t) are connected by the mass conservation law 



8 



mN{t)E{t) = 1, (10) 
where m stands for the average mass of a strand. It follows from Eqs. (9) and (10) that 

^ ' Stt Jo ^Jol- n{t, d,^) ^ ' 

Let 5* be the average concentration of chains (per unit mass) in a (hypothetical) totally 
disentangled network and M{t) the current number of active entanglements per unit mass. 
Because any tight entanglement doubles the number of chains involved in the knot, the 
concentration of chains in the entangled network reads 

S(t) = H, + 2M(i) = 5o[l + Ki)], (12) 

where So = + 2M(0) is the average number of chains (per unit mass) in an equilibrium 
network (before a test) and the dimensionless function 

.(t) = 2^^w:^'^'°> (13) 

characterizes the ratio of the number of loose entanglements that have transformed into the 
active state within the interval [0,t] to the initial number of chains. Combining Eqs. (11) 
and (12), we obtain 

W{t) = ^^[1 + u{t)f r r . sinMA (14) 

^ ' Stt ^ ^ Jo ^ Jo l-n{t,'d,(p) ^ ' 

The average mechanical energy per chain is given by 

. N Wit) 

Substitution of Eqs. (12) and (14) into this equality yields 

^ ' Stt ^ ^ Jo ^ Jo l-n{t,'d,(p) ^ ' 

Our objective now is to calculate the derivative of the function W with respect to time. 
For this purpose, we suppose that the characteristic rate for transition of bonds from their 
flexed to extended conformation substantially exceeds that for changes in the concentration 
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of active (tight) entanglements. This assumption may be explained by the fact that the 
characteristic length of a chain (the length between neighboring entanglements) dramatically 
exceeds that for a strand (the length-scale associated with a bond). This implies that the 
function i/(i) may be treated as a constant when expression (14) is differentiated with respect 
to t. Bearing this hypothesis in mind, we find from Eq. (14) that 



dt 



— n{t, d, Lp) dt 



It follows from Eq. (8) that 



[l-n{t,^,if)Y dt 



(t,T9,ip) sin« + Ji(i), 
(^,1?, ip) sin«. 



(16) 



de, y. 1 
—(t,l) = ^ 

9r ' ^ 2[/-c(t)-/] 



7] 



dn 



1 + Tjnit, I) dt 
Differentiation of Eq. (6) with respect to time yields 
dC 



(tj). 



(17) 



dt 



{t) = Vov{t) ■ [Vof{t)] + Vof(t) • [Vov{t)y, 



where v{t) = dr{t)/dt is the velocity vector for the network. Taking into account that 



Vov{t) = Vof(t) • V{t)v{t), 



where V(t) is the gradient operator in the deformed state at time t, we obtain 
dC 



dt 
where 



(t) = 2Vof(t) ■ D{t) ■ [Vof(t)] ' , 



(18) 



D{t) = - V{t)v{t) + (Wmt)) 



is the rate-of-strain tensor for the network. It follows from Eq. (18) that 

[• —it) ■ r = 2[- Vof(i) • D{t) ■ [Vor{t)r ■ 2F{t,J) : D{t), 

where 



(19) 



F{t,l) ^[Vor{t)V ■{l^l)-^of{t) 



(20) 
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is the generalized Finger tensoiO, the colon stands for convolution and ® denotes tensor 

product. Substitution of Eqs. (17) and (19) into Eq. (16) implies that 

dW 

— it) = T(t) : D{t) - Jit), 

= —A [1 + / dip / — — - . sim^di^, 



47r ' ''Wo ^ Jo 1 -n{t,i^,ip)l .C{t) -l 



4. Constitutive equations 

Observations evidence that under cyclic loading with small frequency (less than 10 Hz), 
the temperature increment is negligible and temperature T remains close to its reference 
value To§. This means that the effect of temperature on material parameters, as well as 
thermal expansion of the network may be disregarded. It is assumed that the deformation 
process is rather slow, which implies that at any instant t > thermodynamic potentials are 
correctly defined. For affine deformation of an incompressible network, the Clausius-Duhem 
inequality reads0 
dQ dT d^ 1^^ ^ 1 - . 



dt dt dt T 

where p is mass density, q is the heat flux vector, is the deviatoric component of the 
Cauchy stress tensor a, \1/ is the free (Helmholtz) energy, S is the entropy and Q is the 
entropy production per unit mass. We accept the following expression for the free energy: 

^ = ^0 + (c - S'o)(T - To) -cT\n^ + W, (23) 

-^o 

where 5*0 and \l/o are the entropy and the free energy in the equilibrium stress-free state at 
the reference temperature To and c is the specific heat. The second and third terms on the 
right-hand side of Eq. (23) characterize the energy of thermal motion. Unlike conventional 
theories of rubber elasticityisi, the terms associated with configurational entropy of chains 
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are disregarded in Eq. (23). According to Ref-Ell, this approximation is acceptable, provided 
that 

H > kBTN, 

where ks is Boltzmann's constant. In the sequel, we suppose that this inequality is satisfied. 

Substituting Eqs. (14) and (23) into Eq. (22) and assuming changes in the concentration 
of junctions to be rather slow (which means that Eq. (21) for the derivative of the mechanical 
energy W with respect to time may be used), we arrive at the formula 

dt \p J \ Tq) dt pT^ - ^ ^ 



Applying the conventional reasoningpl to Eq. (24), we find that the expressions in braces 
vanish. This assertion results in the standard formula for the entropy 

S = So + cln^, (25) 

where the configurational entropy of the network is neglected, and, together with Eqs. (20) 
and (21), the constitutive equation 

a(t) = -P(t)/ + G[l + Kt)f[Vof(t)r- rdv r sinM^ • Vof(t), 

Jo Jo 1 - n[t,'tf,ipj I . C{t) ■ I 

(26) 

where P{t) is pressure and G = pfimE g/(47r). We substitute Eqs. (25) and (26) into Eq. 
(24) and find that the rate of entropy production is nonnegative for an arbitrary loading 
program, provided that 

1. the heat fiux vector q obeys the Fourier law q = — kVT with a positive thermal 
diffusivity k, 

2. the function n{t,-&,(f) satisfies the kinetic equation 

^(t, ^, = a{t)e{t, ^, - <t, ^,^)], n{0, ^, ^) = 0, (27) 

where a is a nonnegative function of time. 
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We postulate that the rate of trans-cis transformation is proportional to the average number 
of strands in a chain, 



where q;i > is a material constant. Equation (28) is explained by the fact that the force 
opposing transition of a bond from its flexed to extended conformation is driven by the 
action of neighboring bonds, whereas stresses in these bonds are inversely proportional to 
the number of strands in a chain, see Eq. (3). Combining Eqs. (27) and (28) and using Eqs. 
(10) and (12), we arrive at the nonlinear differential equation 



with CKo — Q!i/(?7i5o). 

5. Evolution of the concentration of junctions 

To close the model, the kinetics of stress-induced changes in the concentration of temporary 
junctions should be described. Denote by Mq the total concentration (per unit mass) of 
entanglements (both loose and tight). Adopting the first-order kinetics for the mechanically- 
induced evolution of the number of tight entanglements, we assume that the rate of increase 
in the concentration of active entanglements is proportional to the current number of loose 
ones. 



a{t) = aiN{t), 



(28) 





(29) 



dM 



(t)^ f34t)[Mo-M{t)], 



(30) 



dt 



where /3*(t) is a positive function. It follows from Eqs. (12) and (30) that 



dp 
'dt 



{t) = /3.(t)h-Kt)], 



(31) 



where the quantity 



I/O = 



^ Mo-M(O) 
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equals the ratio of the initial number of loose entanglements to the initial number of chains. 
We postulate that the rate of growth for the number of active entanglements, jS^,, is pro- 
portional to the mechanical energy per chain and inversely proportional to the number of 
strands in a chain, 

m = Po^, (32) 

where /3o > is a material constant. The proportionality of /3* to the average mechanical 
energy, wq, reflects the fact that tightening of loose entanglements is driven by mechanical 
factors. In general, any measure of straining may be chosen instead of wq on the right- 
hand side of Eq. (31), see a discussion of this issue inllUll. We introduce the strain energy 
density, wq, by analogy with conventional models for damage of elastomersS. The inverse 
proportionality of /9* to the average number of strands in a chain, A^, reflects slowing down 
of the process of tightening entanglements with an increase in the average length of chains. 
The latter is driven by the growth of the chains' mobility (estimated in terms of the number 
of available configurations). 

It follows from Eqs. (10), (12), (31) and (32) that the function z/(t) obeys the kinetic 
equation 

^(t) = - i^{t)] r r ^ ^'(^'/'f) smMd, z/(0) = 0, (33) 

dV ' 47r^ ^ Jq ^ Jo l-n(t,i9,(/?) > , \ j 

where P = /3oyu/2. Governing equations (26), (29) and (33) are determined by 5 adjustable 
parameters: an analog of the shear modulus G, the rate of trans-cis transition ao, the rate 
of stress-induced increase in the concentration of active entanglements (3, the constant rj 
which characterizes the end-to-end elongation of a chain driven by transformation of a bond 
from its flexed to extended conformation, and the ratio, uq, of the initial number of loose 
entanglements to the initial number of chains. To determine these quantities, we analyze 
uniaxial extension of a rod. 
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6. Uniaxial tension of a specimen 



At uniaxial tension of an incompressible medium, Cartesian coordinates in the deformed 
state, Xi, are expressed in terms of the Cartesian coordinates in the stress- free state, X^, by 
the formulas 

xi = k{t)Xi, X2 = k~^t)X2, X2, = k~^{t)X2„ 

where k = k{t) is the extension ratio. It follows from these equalities and Eq. (6) that 



Vor(t) = k{t)eiei + k 2 (t)(e2e2 + 6363), 
C{t) = A;'(i)eiei + A;-' (i) (6362 + 6363). 

The unit vector I is given by 

I — cos'i^ei + sin'i?(cos(^e2 + sin (^63), 

which implies that the tensor I ®l is determined by the matrix 

cos^ 'd sin 'd cos 'd cos ip sin 'd cos 'd sin (p 
I ®l = sin 7? cos 'd cos </? sin^ 'd cos^ ip siv? ■& sin (p cos (p 
sin 1? cos ■»? sin (p sin^ 1? sin </? cos (p sin^ ■»? sin^ (/? 

Equations (34) and (35) yield 
l.C{t) -1 = k\t) cos^^ + k-\t)sm^^. 



(34) 



(35) 



(36) 



(37) 



Because of the axial symmetry of deformation, the functions n = n{t^'d) and e = e{t^d) are 
independent of which implies that all terms but I ®l on the right-hand side of Eq. (26) 
are independent of </?. It follows from Eq. (36) that 

1 /-Stt ^ _ \ 

— / I ® ld(p — cos^ 'j^eiCi -|- - sin^ '&{e2e2 + 6363). 
2% Jo 2 

Substitution of this expression and Eqs. (34) and (35) into Eq. (26) imphes that 
a{t) = a{t)eiei (7o(i)(e2e2 6363), 
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where 



a{t) = -P{t) + 2ttG[1 + iy{t)f - ^ ^ 

Jo k 



cos2t9 + A;-i(t)sin2'(?l -n(t,i?) 



sin 

sini?(ii?. (38) 



The boundary condition on the lateral surface of the specimen implies that cro(t) = 0. 
Combining this equality with Eq. (38), we find the longitudinal stress 

2k\t) cos^ ^ - k-\t) sin^ ^ e{t,^) 



a{t)^T:G[l + v{t)f r 
Jo 



sin ■&d'd. 



lo k^{t)cos'^^ + k-\t)sm'^^ l-n{t,^) 
Introducing the notation 

z — cos'j?, e{t,z) — e{t,-&), h{t,z) — n{t,'&) 

and bearing in mind that e and n are even functions oi z, we obtain 

2e{t)z^ - k-\t){l - z^) e{t,z) 



a{t) = E[l + u{t)r f 
Jo 



k^{t)z^ + k-\t){l-z^) l-h{t,z) 



dz, 



(39) 



where 
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E = ^pfimEl (40) 
It follows from Eqs. (8), (34) and (35) that 

In the sequel, we focus on stretching with a constant rate of engineering strain, eo > 0, 
k{t) = 1 + eot. (42) 
Combining Eqs. (29) and (42), we arrive at the kinetic equation 



dn 
dk 



277(1 - n) 

1 + 7)71 



— e 



e, n(l, z) = 



(43) 



where a = a/io. Substituting expression (42) into Eq. (33) and calculating the integral over 
(f, we find that 



16 



% = b{uo - ^) f f4, Kl) = 0, (44) 
dk Jo I — n 

where b = P/eQ. Equations (39), (41), (43) and (44) are determined by 5 adjustable pa- 
rameters: E, a, b, 1] and uq. This number is quite comparable with the number of ex- 
perimental constants employed in conventional stress-strain relations in finite elasticity of 
elastomergHlcl Si. 



7. Comparison with experimental data 

To determine adjustable parameters, we fit observations for natural rubber vulcanizates 
in uniaxial tensile tests with the strain rate eo = 2.0 min~^ at room temperature. For a 
description of specimens and the experimental procedure, see Refs.00. 

We begin by matching experimental data in cyclic tests with the maximal elongation 
^max = 6.0 for a vulcanizate with an unspecified composition. First, we approximate the 
stress-strain curve for a virgin sample. Given a Young's modulus, E, the constants a, b, r) 
and uq are found by the steepest-descent algorithm. The parameter E is determined by the 
least-squares technique. Afterwards, we fix the values b and (which are responsible for 
the evolution of the concentration of junctions) and match stress-strain curves measured 
after retraction by using only three constants: E, a and t]. Experimental data in tests with 
various numbers of cycles, i, are depicted in Figure 1 together with results of numerical 
simulation. Material constants are listed in Table 1. The parameters E, a and 1] are plotted 
versus the number of cycles i in Figures 2 and 3, which show that observations are fairly 
well approximated by the stretched exponential functions 

E = Eo-10'^'', a = ao-10'^\ r] = r]o ■ 10'''\ (45) 

where the constants Eq, oq, 170 and ke, i^a, i^q are found by the least-squares algorithm. 

We proceed by fitting observations for a natural rubber vulcanizate with an unspeci- 
fied composition in tensile loading-unloading tests with an increasing maximal amplitude 
of stretching, fcmax- First, we approximate experimental data for a virgin specimen and de- 
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termine adjustable parameters using the steepest-descent procedure. Afterwards, we fix tlie 
constants b and uq and match stress-strain curves measured after retraction with the help of 
three adjustable parameters: E, a and rj. Figure 4 demonstrates good agreement between 
experimental data and results of numerical simulation with adjustable parameters collected 
in Table 2. The quantities E, a and rj are plotted in Figures 5 to 7 versus the number of 
cycles i. These figures reveal an acceptable quality of fitting experimental data by Eq. (45). 

To evaluate the influence of annealing (24 h at the temperature T = 100 °C) on ad- 
justable parameters in the constitutive equations, we approximate observations in tensile 
tests for three vulcanizates: A (polysulfide crosslinks), B (monosulfide crosslinks) and C 
(carbon-carbon crosslinks) with various amounts of crosslinkers (sulfur and dicumyl perox- 
ide, respectively). A detailed description of the composition of specimens in presented in0. 
For any type of vulcanizates, we begin by fitting the stress-strain curve for a virgin sample 
with the minimum content of crosslinker, 0, and find experimental constants by using the 
steepest-descent algorithm. Afterwards, we fix the quantities b and z/q and match other ob- 
servations by using only 3 adjustable parameters: E, a and rj. Experimental data together 
with results of numerical simulation are depicted in Figures 8 to 10 for vulcanizate A, in 
Figures 11 to 13 for vulcanizate B and in Figures 14 and 15 for vulcanizate C. Experimental 
constants are collected in Tables 3 to 5. 

To compare the effects of thermal recovery and recovery by swelling, we match experi- 
mental data in tensile tests for a virgin specimen and for the same specimen after retraction 
(vulcanizates A, B and C with unspecified compositions), and repeat the same procedure 
for a recovered sample (swolen in benzene for 24 h and dried in vacuum). Adjustable pa- 
rameters are found by using the same numerical procedure as for the specimens annealed at 
an elevated temperature. Observations and results of numerical simulation are depicted in 
Figures 16 to 18, whereas material constants are listed in Table 6. 
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8. Discussion 

Figures 1, 4, 8 to 18 demonstrate fair agreement between experimental data and results 
of numerical analysis, which implies that the model may be applied to fit observations in 
uniaxial tensile tests with the axial elongation up to kmax — 8.0. 

Figures 3 and 7 show that the parameter rj increases with the number of cycles, i. This 
conclusion is also confirmed by the data hsted in Tables 3 to 6 (results for the first and 
second stretching). On the contrary, recovery of specimens (by annealing at an elevated 
temperature and by swelling) reduces rj. This means that the quantity rj may be treated 
as a parameter responsible for material damage at the micro-level. Because rj reflects the 
average end-to-end elongation of a chain driven by an individual trans-cis transition (which 
implies that it may be thought of as an average measure of coiling of a chain), an increase 
in rj under cyclic stretching is tantamount to an increase in the average size of a globule 
formed by a chain. This picture results in the conclusion that recovery of natural rubber 
vulcanizates may be treated as coihng of long chains. Tables 3 to 5 evidence that this process 
occurs more intensively for rubbery polymers with high concentrations of crosslinkers, and 
vulcanizates with dicumyl peroxide are annealed more effectively than those with sulfide 
crosshnks. These tables also show that the parameter rj monotonically decreases with the 
growth of the content of crosslinker, (p. Because the growth of the concentration of chemical 
crosslinks results in an increase in the volume fraction of junctions (which is confirmed by the 
growth of Young's modulus with 0), and, as a consequence, a decrease in the average length 
of a chain, N, we may conclude that short chains are coiled up relatively more strongly [the 
qualification "relatively" is important here, because r] is multiplied by the number of strands 
in a chain, N, in Eq. (7)] than longer ones. 

Tables 3 to 6 show that at the second stretching the rate of transition from the flexed 
to extended conformation of bonds, a, exceeds that for a virgin specimen (the data for 
vulcanizate B in Table 6 are the only exception from this rule) . This observation may also 
be treated as some kind of material damage: if we take a bond as a hinge, then straightening 
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two rods linked by the hinge occurs more easily for a preloaded sample than for a virgin 
one. Tables 3 to 5 reveal that a substantially grows with the concentration of crosslinker, 
(j). This increase may be associated with appropriate changes in r]: an increase in the degree 
of coiling makes trans-cis transitions easier. In most cases, annealing of specimens (thermal 
and by swelling) results in a decrease in a which is accompanied by a decrease in 77. However, 
for vulcanizate A with high concentrations of crosslinker, an inverse tendency is observed: 
despite a decrease in 77, the rate of transition from the flexed to extended conformation 
increases, see Tables 3 and 6. 

Young's modulus E strongly increases with the content of crosslinker. This conclusion is 
in accord with Eq. (40) which implies that E is proportional to the square of the initial num- 
ber of chains. Tables 3 to 6 demonstrate that the elastic modulus for a recovered specimen 
exceeds that for a prestrained material. The growth of E during recovery seems quite nat- 
ural, provided that the recovery process is associated with an increase in the concentration 
of entanglements. 

Unlike experimental data collected in Table 2, Table 1 demonstrates that cyclic loading 
results in a pronounced decrease in a and E with the number of cycles, i. Figure 2 reveals 
that this decrease may be fairly well approximated by the stretched exponential law. To 
explain the difference between the data listed in Tables 1 and 2, we refer to Figures 5 to 7, 
which evidence that E and a decrease with i at small amplitudes of stretching and strongly 
increase at large amplitudes. Changes in adjustable parameters observed in the cyclic test 
with an increasing amplitude are similar to those presented in Table 1 during the first 4 
cycles and are analogous to those presented in Tables 3 to 5 when the maximum amplitude 
of stretching, fcmax, exceeds 4. One can only speculate about a physical mechanism for 
this phenomenon. As a possible explanation, we suppose that cyclic loading affects the 
internal structure of a rubber vulcanizate in different ways at small and large amplitudes 
of stretching. Periodic loading with amplitudes which do not exceed some threshold results 
in material training similar to that occurring at isothermal annealing (Table 1). On the 
contrary, cyclic stretching with amplitudes exceeding the threshold strain results in damage 
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of the internal structure which is reflected by the data collected in Tables 2 to 5. 

For each series of tests, the parameters b and vq are found from the condition of the best 
fit for one stress-strain curve and remain fixed in matching observations for other curves. 
Excellent agreement between results of numerical simulation and experimental data leads 
to the conclusion that these parameters are not affected by mechanical factors. 

Table 6 reveals that the differences between the parameters a and E found by fitting 
observations on recovered specimens at the first and second stretching are rather small 
compared with those for virgin samples. This means that recovery by swelling is equivalent 
to equilibration of rubbery polymers. On the contrary, analogous differences in t] remain 
practically constant, which confirms the main hypothesis of this study that stress-softening 
of elastomers may be associated with uncoiling of long chains. 

9. Conclusions 

Constitutive equations have been derived for the isothermal mechanical response of elas- 
tomers at finite strains. Stress-strain relations are developed using the laws of thermody- 
namics and are applied to fit experimental data in uniaxial tensile tests for several natural 
rubber vulcanizates. Fair agreement is demonstrated between experimental data and results 
of numerical simulation. 

We study the effects of cyclic loading, annealing at elevated temperature and recovery 
by swelling on material constants. The following conclusions are drawn: 

1. Stress-softening of elastomers is reflected by changes in the parameter r) which charac- 
terizes a mechanically-induced increase in the average size of globules formed by long 
chains. The quantity rj is increased under periodic stretching with large amplitudes 
(uncoiling of chains) and is decreased at recovery (by annealing or swelling). 

2. In contrast to conventional models for stress-softening of rubbery polymers, fitting of 
observations demonstrates that Young's modulus is altered relatively weakly. This 
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implies that damage should be associated rather with changes in the topology of a 
temporary network than with rupture of long chains. 

3. Cyclic straining results in a noticeable increase in the rate of trans-cis transition of 
bonds, a, with the number of cycles, i (which is fairly well described by the stretched 
exponential law). On the contrary, recovery of specimens (by thermal annealing and 
swelling) induces a decrease in a. 

4. Young's modulus, E, and the rate of transformation of bonds from their flexed to 
extended conformation, a, substantially increase with the content of crosslinker, (f), 
whereas the parameter 77 is a decreasing function of (p. 

5. Adjustable parameters change in a similar way during thermal annealing of specimens 
at an elevated temperature and during recovery by swelling. 

6. Cyclic loadings with relatively small and large amplitudes affect the internal structure 
of elastomers in different ways. Periodic straining with large amplitudes causes dam- 
age to the internal structure, whereas the influence of that with small amplitudes is 
tantamount to recovery of vulcanizate and may be thought of as a kind of training of 
rubbery polymers. 

7. The material parameters b and i/q which reflect the kinetics of mechanically-induced 
transition of loose entanglements into the active state are rather robust. They weakly 
depend on the intensity of loading, as well as on the composition of vulcanizates and 
on the content of crosslinkers. 

8. Recovery by swelling reduces the increments of Young's modulus and the rate of trans- 
cis transition driven by a cycle of loading. This implies that swelling of a rubber 
vulcanizate induces equilibration of its mechanical properties. 

The above assertions are obtained in the framework of the phenomenological model, which 
implies that they should be treated with caution. These conclusions may, however, be 
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taken as rather plausible assumptions (mechanically-induced uncoiling of chains, training of 
an elastomer by sub-threshold periodic stretching, acceleration of trans-cis transformations 
under loading, disentanglement at retraction, etc.) which have to be checked experimentally. 
Vahdation of these hypotheses will be the subject of a subsequent study. 
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Figure legends 



Figure 1: The true stress a MPa versus the engineering strain e for natural rubber with 
an unspecified composition. Circles: experimental datalll. Solid lines: results of numerical 
simulation. Curve 1: i = 1; curve 2: i = 2; curve 3: i = 3; curve 4: i = 8 

Figure 2: Young's modulus E MPa (unfilled circles) and the dimensionless parameter a 
(filled circles) versus the number of cycles i. Symbols: treatment of observationJll. Solid 
lines: approximation of the experimental data by Eq. (45). Curve 1: Eq = 12.3705, 

= -0.0716; curve 2: Oq = 9.2561, = -0.1219 

Figure 3: The dimensionless parameter t] versus the number of cycles i. Cirles: treatment 
of observationslli. Solid lines: approximation of the experimental data by Eq. (45) with 
rio = 4.1975 and = -0.2055 

Figure 4: The true stress a MPa versus the engineering strain e for natural rubber 
vulcanizate with an unknown composition. Circles: experimental dataS. Solid lines: results 
of numerical simulation. Curve 1: a virgin specimen; curve 2: the 4th stretching; curve 3: 
the 5th stretching; curve 4: the 6th stretching; curve 5: the 7th stretching; curve 6: the 8th 
stretching; curve 7: the 9th stretching 

Figure 5: Young's modulus E MPa versus the number of cycles i. Circles: treatment of 
observationsii. Solid lines: approximation of the experimental data by Eq. (45). Curve 1: 
Eo = 21.5638, ke = -0.0125; curve 2: Eq = 3.3807, ke = 0.5003 

Figure 6: The dimensionless parameter a versus the number of cycles i. Circles: treat- 
ment of observations^. Solid lines: approximation of the experimental data by Eq. (45). 
Curve 1: = 9.9095, k„ = 0.0863; curve 2: = 3.7108, k„ = 0.3952 

Figure 7: The dimensionless parameter t] versus the number of cycles i. Circles: treat- 
ment of observations^. Solid lines: approximation of the experimental data by Eq. (45). 
Curve 1: r]o = 4.6700, = -0.0035; curve 2: r]o = 2.7620, = 0.2524 

Figure 8: The true stress a MPa versus the engineering strain e for natural rubber 
vulcanizate A ((f) = 1.25 phr). Circles: experimental data0. Solid lines: results of numerical 
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simulation. Curve 1: the first stretching; curve 2: the second stretching; curve 3: stretching 
after anneahng 

Figure 9: The true stress a MPa versus the engineering strain e for natural rubber 
vulcanizate A (0 = 2.50 phr). Circles: experimental data0. Solid lines: results of numerical 
simulation. Curve 1: the first stretching; curve 2: the second stretching; curve 3: stretching 
after annealing 

Figure 10: The true stress a MPa versus the engineering strain e for natural rubber 
vulcanizate A (0 = 4.17 phr). Circles: experimental dataS. Solid lines: results of numerical 
simulation. Curve 1: the first stretching; curve 2: the second stretching; curve 3: stretching 
after annealing 

Figure 11: The true stress a MPa versus the engineering strain e for natural rubber 
vulcanizate B (0 = 0.2 phr). Circles: experimental data@. Solid lines: results of numerical 
simulation. Curve 1: the first stretching; curve 2: the second stretching; curve 3: stretching 
after annealing 

Figure 12: The true stress cr MPa versus the engineering strain e for natural rubber 
vulcanizate B (0 = 0.4 phr). Circles: experimental data@. Solid lines: results of numerical 
simulation. Curve 1: the first stretching; curve 2: the second stretching; curve 3: stretching 
after annealing 

Figure 13: The true stress a MPa versus the engineering strain e for natural rubber 
vulcanizate B (0 = 0.6 phr). Circles: experimental dataS. Solid lines: results of numerical 
simulation. Curve 1: the first stretching; curve 2: the second stretching; curve 3: stretching 
after annealing 

Figure 14: The true stress a MPa versus the engineering strain e for natural rubber 
vulcanizate C (0 = 0.5 phr). Circles: experimental data@. Solid lines: results of numerical 
simulation. Curve 1: the first stretching; curve 2: the second stretching; curve 3: stretching 
after annealing 

Figure 15: The true stress cr MPa versus the engineering strain e for natural rubber 
vulcanizate C (0 = 3.5 phr). Circles: experimental dataS. Solid lines: results of numerical 
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simulation. Curve 1: the first stretching; curve 2: the second stretching; curve 3: stretching 
after anneahng 

Figure 16: The true stress a MPa versus the engineering strain e for natural rubber 
vulcanizate A with an unspecified composition. Circles: experimental data0. Solid lines: 
results of numerical simulation. Curve 1: a virgin sample, the first stretching; curve 2: 
a virgin sample, the second stretching; curve 3: a recovered sample, the first stretching; 
curve 4: a recovered sample, the second stretching 

Figure 17: The true stress a MPa versus the engineering strain e for natural rubber 
vulcanizate B with an unspecified composition. Circles: experimental dataS. Solid lines: 
results of numerical simulation. Curve 1: a virgin sample, the first stretching; curve 2: 
a virgin sample, the second stretching; curve 3: a recovered sample, the first stretching; 
curve 4: a recovered sample, the second stretching 

Figure 18: The true stress a MPa versus the engineering strain e for natural rubber 
vulcanizate C with an unspecified composition. Circles: experimental data0. Solid lines: 
results of numerical simulation. Curve 1: a virgin sample, the first stretching; curve 2: 
a virgin sample, the second stretching; curve 3: a recovered sample, the first stretching; 
curve 4: a recovered sample, the second stretching 
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TABLES 



Table 1. Adjustable parameters for NR vulcanizate 



a r] E MPa uq b 

1st stretching 8.9 4.12 12.48 1.7 0.6 

2nd stretching 8.0 4.41 10.89 1.7 0.6 

3rd stretching 7.8 4.50 10.13 1.7 0.6 

8th stretching 5.5 4.57 8.79 1.7 0.6 



Table 2. Adjustable parameters for NR vulcanizate 



a r] E MPa uq b 

1st stretching 9.9 4.67 21.65 1.7 0.35 

4th stretching 7.9 4.36 14.69 1.7 0.35 

5th stretching 14.1 4.63 19.52 1.7 0.35 

6th stretching 14.4 4.93 20.85 1.7 0.35 

7th stretching 17.0 5.19 23.61 1.7 0.35 

8th stretching 19.2 5.53 29.64 1.7 0.35 

9th stretching 23.2 5.95 41.66 1.7 0.35 
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Table 3. Adjustable parameters for vulcanizate A 





a 


V 


E MPa 


^0 


b 


(f> = 1.25 plir 


1st stretching 


13.5 


4.27 


12.66 


1.6 


0.4 


2nd stretching 


15.1 


4.67 


10.25 


1.6 


0.4 


after annealing 


15.0 


4.47 


11.41 


1.6 


0.4 


= 2.50 phr 


1st stretching 


11.5 


3.55 


15.15 


1.6 


0.4 


2nd stretching 


14.4 


4.11 


15.15 


1.6 


0.4 


after annealing 


18.2 


3.90 


17.78 


1.6 


0.4 


= 4.17 phr 


1st stretching 


12.3 


3.21 


22.24 


1.6 


0.4 


2nd stretching 


17.6 


3.32 


18.66 


1.6 


0.4 


after annealing 


18.4 


3.32 


21.06 


1.6 


0.4 
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Table 4. Adjustable parameters for vulcanizate B 





a 


V 


E MPa 


^0 


b 


= 0.2 phr 


1st stretching 


10.9 


4.60 


13.43 


1.7 


0.2 


2nd stretching 


16.8 


5.10 


14.35 


1.7 


0.2 


after annealing 


16.4 


4.81 


16.04 


1.7 


0.2 


(j) = 0.4 phr 


1st stretching 


14.4 


3.90 


22.69 


1.7 


0.2 


2nd stretching 


24.7 


4.20 


24.26 


1.7 


0.2 


after annealing 


20.4 


4.00 


25.01 


1.7 


0.2 


= 0.6 phr 


1st stretching 


20.4 


3.66 


34.84 


1.7 


0.2 


2nd stretching 


34.7 


3.83 


35.40 


1.7 


0.2 


after annealing 


33.0 


3.62 


39.12 


1.7 


0.2 



31 



Table 5. Adjustable parameters for vulcanizate C 





a 


V 


E MPa 


^0 


b 


= 0.5 phr 


1st stretching 


8.4 


5.61 


8.37 


1.7 


0.3 


2nd stretching 


11.7 


6.02 


7.01 


1.7 


0.3 


after annealing 


11.1 


5.65 


8.65 


1.7 


0.3 


= 3.5 phr 


1st stretching 


28.2 


3.58 


23.33 


1.7 


0.3 


2nd stretching 


38.1 


3.67 


22.34 


1.7 


0.3 


after annoalino 


35.9 


3.57 


23.82 


1.7 


0.3 
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Table 6. Adjustable parameters for natural rubber vulcanizates 





a 




E MPa 




b 


Vulcanizate A 


virgin samples 


1st stretching 


10.9 


4.40 


23.48 


1.6 


0.4 


2nd stretching 


16.1 


5.11 


24.29 


1.6 


0.4 


recovered samples 


1st stretching 


16.7 


5.01 


24.11 


1.6 


0.4 


2nd stretching 


16.1 


5.25 


24.38 


1.6 


0.4 


Vulcanizate B 


virgin samples 


1st stretching 


12.5 


5.64 


31.97 


1.6 


0.2 


2nd stretching 


9.5 


6.04 


23.52 


1.6 


0.2 


recovered samples 


1st stretching 


8.6 


5.63 


24.21 


1.6 


0.2 


2nd stretching 


10.5 


5.91 


23.62 


1.6 


0.2 


Vulcanizate C 


virgin samples 


1st stretching 


7.8 


5.62 


25.34 


1.6 


0.2 


2nd stretching 


11.6 


6.07 


26.47 


1.6 


0.2 


recovered samples 


1st stretching 


9.6 


5.66 


27.21 


1.6 


0.2 


2nd stretching 


9.6 


6.23 


25.55 


1.6 


0.2 
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